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SUMMARY 

By using Galerkin’s variational method we build up an approximate solution 
for a system of two differential equations with linear partial derivatives of 
the second order. This system of differential equations corresponds to the 
physical model, known in the literature as the Timoshenko Beam. The results 
obtained can be finally applied to two particular cases representing respectively: 
the case of a beam with a rectangular section, with a constant height and a basis 
with a linear variation: 

the case of a beam with a constant basis and a height with cubic variation. 

INTRODUCTION 

We are taking into consideration a heterogenous elastic straight beam pos- 
sessing variable geometrical and mechanical characteristics all along the beam. 

We are considering the small, cross-cut non-damping forced oscillations. 

The mathematical model chosen to be subjected to analysis consists in a system 
of two linear equations with partial derivatives of second order, corresponding 
to the physical model known in the literature under the name of Timoshenko Beam. 
This model is more exact than the classical one usually employed in the engineer- 
ing calculations, that is the Euler-Bernoulli model. The difference between them 
consists in the fact that while for the Euler-Bernoulli model only the deforma- 
tions given by the bending moment or by the translation inertia are taken into 
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account, in the Timoshenko model the transverse shear and the rotational inertia 
are also taken into consideration. As a result the Timoshenko model reflects 
more exactly the physical reality. It is well-known that (ref. 1) the differ- 
ences between the two theories become significant in the case of (relatively) 
short beams and this cannot be neglected any longer. 

Although the literature referring to the dynamics of the Timoshenko Beam 
is abundant enough, the matter of the non-damping beam has been insufficiently 
treated. 

In the present paper we try to determine the approximate solutions of the 
phenomenon by means of the Galerkin variational method. We are of the opinion 
that the above mentioned method is most suitable in solving the subject consid- 
ered. The choosing of the system of coordinates required by the Galerkin method 
assures the convergence of the obtained solutions. 

SYMBOLS 

6 (x-^ ) Dirac function 

/X 

x^ ^(l-x)® dx 
u 

K coefficient of the form of the section 

G cross-cut modulus of elasticity 

p density of material 

E longitudinal modulus of elasticity (Young) 

A(x) area of cross-cut section 
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I(x) 


moment of inertia of cross-cut section 


W(x,t) 


cross-cut displacement 


4^(x,t) 


rotation angle 


(f(x),g(x)) scalar product: (f,g) = J' f(x)g(x)dx 


I 

V(x) 

U(x) 

a,X 


c;[o:,i] 


length of beam 


time-independent cross-cut displacement 


time- Independent rotation angle 


cross-sectional area parameters 


moment of inertia parameter 


class of functions defined on 0 to 1 


THE DIFFERENTIAL EQUATIONS OF THE FHENOWENON 
The differential equations for the phenomenon are as follows: (ref. 2) 


ne axrrerenciai equations xor tne 




( 1 ) 


Solutions for the differential equations are determined as follows! 


( 2 ) 
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for boundary conditions 


^ Ci ^ <P (3) 

and f(x,t) is a perturbance force, a mobile, but concentrated force for a unit 
magnitude: 

By considering equation (2) the system of equation (1) becomes two differ- 
ential equations of the fourth order for V(x) and U(x) as follows: 


The differential equations (4) and (5) for the following two cases are as 
follows : 

easel: ^ 1 C*0) ^ 

C5^h<o"^-v cx U ^-v »oX--v 4- 

Q-»5\X4 
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( ‘i A ^ \ / , 

al^xva« ) v + 

^ X /' 

-V 0=^oV-<r^^(VA -V (7) 

4- Q^'cio'^-V ^ <isv^+ '^s-A' 4-V V — 

Case 2: Kc-o\=^ ^Cv-vNx'^ kxA= Kv+>-x^ 

0^s^-v'i?v'v.-vO^{'. Jt •!?<». + 'iy\+0«')L) v0j;i'>i!Afll4^^ A 

<1^XA Css^ U = 

— (SA^J^+Ql^x+'J^^^SV^;VQ*'AA+'^l^'^SAx-^;^ 


C'^XiJX +‘^‘J\X+^iX-V ‘^»?.+'^^+CS<jxi^+(J<5^4 Q-<5V^+ -V 

+ a^k+o<;x^+ ao^+ aaxk "V 

Ak^oX+'k^'+ko^ A" 

4(ii(i^+'o^'x^+Wtti5-v t)t»?^+-'=’<iAXJr k's'xA W<;') x/ -V 
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— V Q + ^y^^^%{l k (X%(^ 

•4“ 





4 . 



+a^^v4%x-vaK\^¥a^\^ 

At v'o ^ Vta-T,^ 


THE APPROXIMATE SOLUTION 

We shall integrate the differential equations (6) , (7) , (8) , and (9) by 
means of the Galerkin method. 

In the case of boundary conditions of equation (3) we shall consider l-unit 
which is always possible by 

X ■• 04^4^-^ 0 4X.<1 

Using the Galerkin method, we shall determine an approximate solution for 
equation (6) as follows: 


I 


(10) 


We choose 4>. (x) of the form (ref. 3) 

iC 


X v5x\ = 'X. (Vt i 
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The system of coordinate functions ^^.s to satisfy the boundary 

conditions of equation (3) which become equivalent with the following 
conditions: 


U \J G t-, <o 


The approximate solution (10) becomes: 



^ \y ^ 


ya-V. 


( 11 ) 


The OL constants are determined out of the following algebraic system: 



( 12 ) 


where L is the left part of equation (6) , and 
same equation, that is: 


g is the right part of the 

r X 



4- V 


) 
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The system of equation (11) is a non-damped algebraic system of n equa- 
tions with n indeterminates. This system is compatible because the determi- 
nant formed with the coefficients of the undeterminants is a Gramm determinant 
of a linear Independent system of functions. For the calculation of the scalar 

product 3 nd we have kept in view the following points; 

K j j 

We have used the Euler's Beta function 

We have used the following fomula (ref. 4) in calculating the scalaf 
product; 

if e 

O 

For equation (7) with the boundary conditions of equation (3) which mean 
V(o)®V(l)==o we are going to give an approximate solution of the following form 

(13) 


/ i> \C 
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where the constant is drawn from the following algebraic system; 

where L is the left side of equation (7) and g is the right side of the 
same equation. 

Analogous to equation (8) we build up an approximate solution of the follow- 
ing form: 

= • 'V(.='^-+s 


where the constants are determined from the following algebraic system: 




where L is the left side of equation (8) and g the right side of the same 
equation. 

Finally, for equation (9) we build up a solution of the following form; 




^ ' ■1- vc"^ C V- X ) > 


where the 6, constants are determined from the following algebraic system: 



V 






where L and g are the left side and right side of equation (9). 
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As a conclusion to case 1 the approximate solutions built up by the Galerkin 


method are the following: 



and, for case 2 the approximate solutions are the following 



O Vt. OAl'U. 

= 0. 'X,0-X^ 


PECULIAR CASES 


In the following lines we shall use the obtained solution for two particulai 
cases, which will be also an Indirect checking of the accuracy of the obtained 
results. 

We build up the first two approximations ^2 respectively W^; 

^2 for following situations: 




They represent respectively the case of a beam with a rectangular section, hav- 
ing a constant height and a base with a linear variation, and the case of a 
beam with a constant base and a height with a cubic variation and this because, 
from an applicative point of view the beam sections are in many cases considered 
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rectangular. Case (a) The equation (6) , if we consider (14) is reduced to the 


following equation 

The first and second approximations are respectively; 


If we compare with 1^2 ^ rectangular beam made of steel we come 

to the conclusion that the two approximations are comparable: = ^2 

certain X values and for certain x values 


X 

0,1 

0,2 

0,3 

0,4 

0,5 

X 

0,252 

0,541 

0,528 

0,573 

0,525 


This conclusion results from the following calculation: 


tX 







I 



The ^i>^2 determined from the following algebraic system:. 
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where 




and for the steel in S.I. units 



Equation (7) 


Q =i ^ X Vg \ t X \(0 ’ Va \ 'trt , 

■ /y^L ) 'VA.t ) ^ 

\o"'Vl. \ ‘ Ko*<S^oV \ CXc^ 

•va’* ’ ^ “ \% 

then becomes: 


K W \ Wi n/ "-V Cn/') =i [_4 , ^\+\xV 4 

-V S Vx-^ \ 

where CX E Q k K ^ “’j W •= C E 4k.cj\ ^ 

e= <^cJ'Q^c*a I^-kQKJ^Kt *) \t 
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The first and second approximations are 



here 




X^= 
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1 


l^va.3i 




"?, _ 1 


1 ^'“j. i - '- 


s - -T- ^'+ ^ >•'+ i-/ ^ ^ 

^ ^ \h\o %K<o 2A I'o \A«o 


^A- 

1 


L 

X^4 



I 



?x4 







-d 

it 

— 

C&A 



>? A 

\ 

-L. X 



A- 

1 L >.4^ 

4-'^ A 

,L 


-( 

■ x> 


^5- 


7j'«^ 


Bec.„. \44.=^i ^U 

is made that the two approximations are comparable for certain X 
for certain x values such as 


the conclusion 
values and 


X 

0,1 

0,2 

0,3 

0,4 

0,5 

X 

0,675 

0,515 

0,585 

0,525 

^ 1 

0,567 


CONCLUSIONS TO THESE PECULIAR CASES 


For equations (8) and (9) we come to the same result, that is: the first 

two approximate solutions are equal for the given values of X for the same 
value of X :0,5: 0 < x 5 1 that is, the approximate solutions are comparable 
among themselves in the vicinity of where the concentrated perturbance force is 
applied: when x = ? = -j* 
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